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ABSTRACT 

Single particle inclusive energy spectra in e^e~ annihilation are analyzed in 
terms of moments. By assuming Local Parton Hadron Duality (LPHD), exper- 
imental data in a wide cm. energy range from 3 GeV up to LEP energy are 
compared to the theoretical predictions of Modified Leading Log Approximation 
(MLLA) of QCD with and without taking into account the running of as ■ MLLA 
with running coupling (Limiting Spectrum) is found to reproduce experimental 
results very well, while the model with fixed coupling is inconsistent with data. 
Rescaled cumulants are shown to be sensitive to the running of as in the asymp- 
totic regime, while the Lorentz- invariant distribution, Edn/(Pp, points out this 
effect at very small energy E of few hundreds MeV. These results give a direct 
evidence of the running of the QCD coupling in inclusive energy spectra and 
lend further support to the LPHD picture. 



1. Introduction 



The main open problem in strong interaction phenomenology is the understanding 
of the hadronization mechanism. There is indeed a gap between theoretical predic- 
tions coming from perturbative QCD at parton level and experimental results at 
final particle level. A bridge between these two sides is given by a soft confinement 
mechanismBi: single particle inclusive spectra are simply required to be proportional 
to single parton inclusive spectra obtained from perturbative QCD, provided the par- 
ton cascade is evolved down to a cutoff scale of the order of the hadronic masses 
Qq ~ nih. This hypothesis of Local Parton Hadron Duality (LPHD)i has been par- 
ticularly successful in describing both the shape and the position of the maximum 
of single inclusive charged particle momentum spectra in hadronic jets in e~^e~ an- 
nihilation ranging from PETRA/PEpi to LEP energyoQ'Qo. Recently, similar results 
have been obtained in DIS for inclusive momentum spectra in the Breit frame at 
HERAi§. 

In this paper, I report on the results obtained in collaboration with Wolfgang 
Ochs0. We performed the analysis of moments up to order 4 of single particle inclusive 
energy spectra in hadronic jets in e~^e~ annihilation. Since moments are sensitive to 
the tail of the distribution, i.e., to the low energy region, a consistent kinematical 
scheme has been defined to relate parton and hadron spectra. Our aim was two-fold: 



1 



firstly, to investigate the sensitivity of single particle energy spectra to tlie running 
of the coupling a^, in particular for soft particles in the low energy region, where 
the coupling is stronger; secondly, to test the validity of the description in terms of 
perturbative QCD plus LPHD in the small cm. energy domain. 

2. The theoretical framework 



2.1. MLLA Evolution Equations 

The evolution equation of single parton inclusive energy distribution in the frame- 
work of MLLA is given byil: 

-^xDf(x,iogP0)= E f'dz^^^^Az)\-DE(-AogizPe)) 

dlOgfc) C^q,g-'^ ^'^ VZ \Z J 

with the boundary condition at threshold P0 = Qo, a;Z)^(x, log Qo) = <^(1 ^ x)5^. 
Here P and denote the primary parton energy and jet opening angle respec- 
tively, X is the energy fraction carried by the produced parton, ^l{z) are the parton 
splitting kernels and A, B, C label quarks, antiquarks and gluons. The QCD 
running coupling is given by its one-loop expression agikt) = 27i/b\og{kt/ A) with 
b = {llNc — 2nf)/3, A the QCD-scale and and nj the number of colors and of 
flavors respectively. The scale of the coupling is given by the transverse momentum 
kt — z{l — z)PQ. The shower evolution is cut off by Qq, such that kt > Qo- 
The integral equation can be solved by Mellin transform: 

D^{Y,X) = ['—x-[xD{x,Y,X)] = r d^e-^^D{^,Y,X) (2) 
JO X Jo 

with Y = log ^ log ^5 A = log X' ^ ~ ¥ parton energy k. 

In flavor space the valence quark and (±) mixtures of sea quarks and gluons evolve 
independently with different "eigenfrequencies" . At high energies, the dominant con- 
tribution to the inclusive spectrum comes from the "plus" -term, which we denote 
by D^{Y, \) = D^{Y,\). By keeping now, as usual, only the leading singularity in 
tu-space plus a constant term, one gets: 

(3) 

where a = llNc/3 + 2nf/3N^ accounts for recoil effects. 
By defining the anomalous dimension according to: 



D^{Y,X)=D^{Yo,X) expj^^' 



dyiu[asiy + \)]] (4) 
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the evolution equation for the inclusive spectrum can also be written as a differential 
equation for 7^. In general this nonlinear equation has two roots and the solution 
consists of a superposition of two terms of the type (Q). The explicit solution of eq. (H) 
can be expressed indeed as a linear superposition of two hypergeometric distributiong3. 
An interesting case of the general solution is the Limiting Spectrum, where the two 
parameters coincide, Qo = A, i.e., A = 0. In this case, the expressions simplify and a 
simpler integral representation of the spectrum has been achieved. 



2.2. MLLA with fixed as 

In order to study the sensitivity of the single particle inclusive energy distribution 
to the running of the coupling, we built a new model based on MLLA formalism but 
with the coupling kept fixed at a given value. In this case, the value of the coupling 
or, equivalently, the value of the parameter 70, 7o = Gas/n, replaces A as a free 
parameter of the model. Accordingly, eq. @ becomes now: 

^D^{Y, A) - ^'MY, A) = -2r^ (^u + D^{Y, A) (5) 

with r] = Q'^q/SNc = aag/'^'n. 

One can then perform explicitly the inverse Mellin transform and get an explicit 
solution for the inclusive energy spectrum in MLLA with fixed-ag: 




3. Moment analysis: from momentum spectra to energy spectra 



3.1. Moment analysis 

In addition to the standard study of single particle inclusive distribution, one 
can also study the multiplicity A/" and the normalized moments of the distribution, 
< e(>",A) >, given by: 



Af = jdiD{i,Y,\) (7) 

<e{Y,X)> = ^Jd^^W{^,Y,X) 

or the cumulant moments Kg{Y, A)Eilii; the latter are related to ordinary moments by 
a standard cluster expansion: 

Ki = <^>=^ , K2=a^=<{^-0'> (8) 
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One also introduces the reduced cumulants kg = Hq/cr'^, in particular the skewness 
s = and the kurtosis k = k^. The cumulants be derived from the expansion: 

logD^iY,X) = J2K,{Y,X)^—^ ^ Kg{Y,X)=\^-—j \ogD^{Y,X) 

withAfE{Y,X)=D^{Y,X)l^,. 

At high energies, one term of the form (^) dominates and one obtains 



, (9) 

LJ = 



(10) 

w=0 



Kg{Y) = Kg{Yo) + J^^ dy y-Q^j 7^[«s(z/)] 

This equation shows the direct dependence of the moments on asiY). For fixed a^, 
for example, one obtains directly i^qiY) oc Y for high energies. 

Let us stress that the study of moments and cumulants is worth to be pursued, 
since explicit analytical expressions for their dependence on cm. energy have been 
obtained^; moments of order q in the MLLA model with running coupling can indeed 
be expressed as a function of a, 6, Y + X and A, while moments for the Limiting 
Spectrum can be expressed in a simpler way in terms of the parameter B = a/b and 



the variable z = ^l^N^Y /h. In the fixed model, the solution is straightforward, 
since in this case the evolution equation for the anomalous dimension reduces to an 
algebraic equation. In addition, if the higher order cumulants (g > 2) are sufficiently 
small, one can reconstruct the ^-distribution from the first four order cumulants, using 
a Distorted Gaussian parametrization of the spectrumlll. 

However, let us notice again that high order moments are sensitive to the tail 
of the distribution; therefore, one has to look carefully at the treatment of the low 
energy region E Qq. This point is discussed in detail in the next Subsection. 

3.2. The resettling procedure 

In theoretical calculations, partons are taken as massless, Ep = pp, but a pt cutoff 
Qo, Pt > Qo, is needed for infrared regularization; therefore, Ep > Qq and the theoret- 
ical spectrum goes linearly to for E ^ Qq. These predictions are usually compared 
to experimental results on single particle inclusive momentum spectrum, which show 
a finite tail down to very low momentum (large ^). The difference of the two spectra 
in their tails could affect significantly the moment analysis. In order to cure this 
anomaly, let us remind that observable hadrons are massive, therefore Eh > mh- It 
has been therefore suggested to identify the cutoff parameter Qo appearing in theo- 
retical calculations with an effective hadron mass averaged over all charged particles. 
In addition, one has to look at energy spectra rather than at momentum spectra. Let 
us also remind that the Lorentz- invariant distribution E-^ approximately follows an 
exponential law in low-energy experiments for sufficiently small energies: 

E^-e-^/^" (11) 
d-^p 
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with Eq ^ 150 MeVlllllllll. By performing a suitable rescaling, the spectrum Ehdn/dph 
goes then linearly to for E Qq. 

In conclusion, we propose to compare theoretical predictions for single particle 
inclusive energy spectra with the experimental spectrum 

vs. f, = log± = log^ (12) 

dph Xe E 

where an effective mass Qq is assigned to all particles, i.e.. Eh = \Jp\ + Ql- Notice 
then that experimental spectra do depend on the chosen value of Qq. 

4. Results 



4^.1. Moments 

Figure 1 shows (the logarithm of) the average multiplicity Me and the first four 
cumulants of charged particle energy spectra as a function of y = log(v^/2(5o) 
for Qq = 270 MeV. For the calculation of moments, an additional value for the 
bin in the unmeasured interval near C,e — Y (small momenta) has been added by 
linear extrapolation. Errors of the moments are taken as the sum of statistical and 
systematic errors. The latter particularly affects the overall normalization and then 
the average multiplicity; higher order moments are on the contrary independent of 
the overall normalization. The errors of the moments also includes the errors on the 
central values of ^e in each bin, taken as half the bin-size. The value of Qo has been 
obtained by comparing the moments determined for a selected Qo with the theoretical 
predictions of MLLA with running at different values of A. The best agreement is 
obtained for the Limiting Spectrum 

Qo = A ^ 270±20 MeV (13) 

Predictions of the Limiting Spectrum with this value of Qo and number of flavors n f 
= 3 are shown in the Figure (solid lines). For the average multiplicity we added two 
new free parameters: 

ATe = ci^2AfLs + C2 (14) 

where the factor 4/9 is the proportionality factor between quark and gluon jets and 
the factor 2 accounts for the two hemispheres; the two additional free parameters 
Ci and C2 give the proportionality factor between parton and hadron spectra and the 
leading particle contribution respectively. This second parameter is needed to recover 
the right behavior near threshold, as the prediction of the Limiting Spectrum does 
not fulfill the correct boundary conditions. The values of the two parameters have 
been fixed by fitting the lowest and the highest energy data points. 
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Fig. 1. The average multiplicity jVb and the first four cumulants of charged particles' energy spectra 
Edn/dp vs. ^£;,i.e., the average value ^b, the dispersion a"^ , the skewness s and the kjttctosis fc, are 
shown as a function of ji" = \og{^/s/2Qo) for Qo = 270 MeV. Data points from MARK iN at ^ = 3, 
4.03, 7.4 GfiV, TASSCH at ^/i = 14, 22,|35, 44 GeV and LEP at Vs = 91.2 GeV (weighted averages 
of ALEPBEI, DELPHIeI, Lai and OPALeI). Predictions of the Limiting Spectrum (i.e. Qo = A) of 
MLLA with running as and of MLLA with fixed as arc also shown (for Uf = 3). Predictions of the 
average multiplicity refer to the two-parameter formula A/b — ci2^Npart + C2. 
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The predictions from MLLA with running and normahzation at threshold are 
remarkably successful considering the fact that there are only two parameters, actually 
coinciding, for the four cumulants. The deviation in at the lower energies may 
suggest that the leading valence quark, which is now neglected by restricting to the 
spectrum, plays indeed a role at low cm. energies. Otherwise, the predictions 
are confirmed at an almost quantitative level down to the lowest cm. energies. 

Figure 1 also shows the theoretical predictions for MLLA with fixed as (dashed 
lines) for the same value of Qo and Uf and 70 = 0.64. This value of 70 has been chosen 
in order to obtain a good fit for the average multiplicity, where we proceed as above 
to a fit with two additional free parameters. With this value of 70, the asymptotic 
slope for is also well reproduced. An adjustement of the absolute normalization 
of ^E at a particular energy Yq would be possible if the perturbative QCD evolution 
towards lower energies and the normalization at threshold are abandoned. Looking 
at higher order moments, however, the model with fixed is unable to reproduce the 
experimental behavior. In addition, the fixed as regime can be excluded already close 
to threshold. Indeed, if we suppose that the coupling is fixed for a certain energy 
interval near threshold Yq < Y < Yi and only runs for Y > Yi, then ^ would be 
shifted towards smaller values because of the very different evolution near threshold, 
in contrast to the experimental behavior. 

4-2. Resettled eumulants 

In addition to the standard moment analysis, it is interesting to consider the 
rescaled cumulants /tg/^; in the high energy limit they become energy independent in 
the model with fixed coupling (see eq. ([T0|)) and therefore exhibit more directly the 
difference to the case of running coupling in the asymptotic domain. 

Figure 2 show the experimental results on rescaled cumulants of order 2, 3 and 4 
for Qq = 270 MeV; a clear F-dependence in the available energy range is visible; this 
result is well reproduced by predictions of MLLA with running and = 3, while 
is again in contradiction with the constant behavior expected for fixed a^. It would 
be interesting to continue this type of studies for instance at TEVATRON, where jets 
of higher energies are available, in order to test whether the rescaled cumulants show 
the remarkable dependence on Y predicted by the theory. 

Let us also comment the dependence of the above results on the number of active 
flavors Hf. In leading order nf enters only through the running coupling as{Y,nf), 
while at next-to-leading level also through the parameter a. It has been argued that 
the number of flavors Uf should be kept fixed at 3, as light quarks should dominate 
quark-pair production in a well-developed cascade. However, one can also include 
heavy quark thresholds in the calculation simply by taking July) under the integral 
with as the number of open flavors at energy y. Restricting our discussion to the 
leading term we note that in the original eq. (|l]) the scale of is kt — z{l—z)P < P/4. 
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Fig. 2. Rescaled cumulants Kg/f as a function of y = log(v^/2Qo) with the corresponding pre- 
dictions of the Limiting Spectrum of MLLA with running as either without or with heavy flavors 
included. Data as in Figure 1. For fixed these quantities approach constant values at high 
energies. 



Therefore, we included the effect of heavy quarks in rij for momentum P > Amq. 
The effect of including heavy quarks thresholds in the model gives a correction of a 
few %, as shown for rescaled cumulants in Figure 2 (dashed lines). In addition, let 
us remind that further uncertainties on the scale of the running, two-loop effects and 
a possible explicit mass dependence of the coupling^ are neglected in this approach. 

4-3. Low- energy behavior: Boltzmann factor or running-ag effect 

Let us consider again the Lorentz-invariant distribution Edn/d^p as a function of 
the particle energy E for small energies up to 1 GeV. The experimental exponential 
behavior is often related to the Boltzmann factor in the thermodynamical description 
of multiparticle production. 

Figure 3a shows experimental data on the invariant distribution at three differ- 
ent cm. energies. All three distributions tend to a common behavior at very low 
energies E, even though deviations from a simple exponential are visible at larger 
cm. energies. Theoretical predictions for the Limiting Spectrum of MLLA with run- 
ning as and for MLLA with fixed as (with 70 = 0.64) with the same value of Qo 
are also shown. Both theoretical curves are normalized to the experimental average 
multiplicity. Figure 3b compares the same data with theoretical predictions from 
the Limiting Spectrum of MLLA with running as in a wide energy range up to 30 
GeV. It is remarkable that the Limiting Spectrum of MLLA is in good agreement 
with experimental data in a wide energy range and can approximately reproduced 
the exponential decrease of the spectrum in the low energy region. On the contrary. 



8 



10 

E— 

[GeV-2 



. 1 



0.01 r 



.001 




0.25 



. 5 



. 75 



. 001 
1 0.25 



0.5 



. 75 



E [GeV] E [GeV] 

Fig. 3. Invariant distribution Edn/cPp as a function of the particle energy E for Qq — 270 Me\L. 
a) comparison of experimental data from MARK W3 (diamonds), TASSOeI (triangles) and OPALQ 
(squares) with predictions from MLLA with running (Limiting Spectrum) (solid lines) and with 
fixed as (dashed lines), b) same as in a), but in the energy range up to 30 GeV; only predictions 
from MLLA with running as are shown; c) analytical predictions of the Born term of DLA with 
running as (solid line) and with fixed as (dashed line); both curves are arbitrarily normalized to a 
common value. 



the exact solution of the fixed-cts model shows a flatter behavior, inconsistent with 
experimental behavior. 

This behavior near the kinematical limit E Qq can be qualitatively understood 
in an analytical way in the DLA, where eq. (|l]) becomes: 

r« . rv^ ,Ca 



D{^, y) = m + / d^' r dy'^^l{y')D{i\ y') 
Jo Jo JMc 



(15) 



where y = log kQ/ Qq = log A;t/ Qq = Y - C.. 

The above evolution equation can be solved iteratively; with one iteration, one 
gets the Born term, which is the leading term in the limit E ^ Qq: 

Ca fy 



(16) 



The difference between the fixed and the running as case is fully due to the expression 
of the anomalous dimension: in the fixed-a^ model, 70 becomes a constant, then: 



DM^,Y) = 6{0 + ^l'oiy-0 



(17) 



whereas for running a^, 7o(l/') = 12/6/(?/' + A) 

n 



(e,l^) = ^(0 + ^ylog^^ 
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Notice that in both cases the spectra go to finite values independent of cm. energy 
in the hmit E ^ Qq; the dependence of the spectrum on cm. energy enters in higher 
orders terms. Predictions for the Born terms of DLA in the two models (with arbitrary 
normalization) are shown in Figure 3c; notice the relative enhancement at low energy 
E in the model with running a^. This behavior can be intuitively explained: as for 
decreasing particle energy E also the typical particle pt is necessarily decreasing, the 
coupling asipt/-^) is rising in the running ag case, whereas the fixed as model cannot 
account for such an effect. 

5. Conclusions 

The analysis of energy spectra in terms of moments has been performed. A consis- 
tent kincmatical scheme to relate parton and hadron spectra in the low-energy region 
has been defined. This yields a good agreement of experimental data with theoretical 
predictions of MLLA with running plus LPHD in a wide cm. energy range from 
3 GeV up to LEP energy. The best agreement with data is found for the Limiting 
Spectrum, where the two parameters of the theory coincide, i.e., Qo = A; the best 
value of Qo has been estimated to be 270±20 MeV. The zero-th order moment, i.e., 
the average multiplicity, requires two more parameters to fix the overall normalization 
for particle production and the leading particle contribution. Predictions of MLLA 
plus LPHD but with fixed coupling arc inconsistent with the experimental behavior 
in the full energy range, thus showing a direct evidence for the running of the cou- 
pling in hadronic energy spectra. The effect of the running as is most pronounced 
near threshold where the variation of the coupling is strongest; it approximately re- 
produces the exponential dependence of the Lorentz-invariant distribution Edn/d^p 
at low particle energies E of few hundreds MeV. The analysis of rescaled cumulants 
is shown to be more sensitive to the difference with the fixed-«s regime in the high- 
energy region, for instance at TEVATRON energies. A similar analysis ranging from 
low to high values of the available energy can be performed also at HERA by looking 
at hadronic jet production in the Breit frame. 

The obtained results give further support to the LPHD picture. A word of caution 
is appropriate nevertheless. The MLLA is based on high energy approximations 
(restriction to contribution, neglect of Next-to-Next-to-Leading effects) which 
have been extended into the low-energy region. A deviation in ^ at 3 GeV could indeed 
suggest that further contributions are needed. Also additional scheme-dependence of 
the coupling and two loop effects have been neglected in the calculation. Finally let 
us remind that we safely identified Qo with an effective hadron mass, since we referred 
to single inclusive spectra for all charged particles, but this identification cannot be 
pushed forward to single inclusive spectra for identified particles without taking into 
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account the effects of resonance decay, as recently shown in 



6. Acknowledgements 

I am grateful to Wolfgang Ochs for fruitful suggestions on the content of this 
note. I would like to thank Ladislav Sandor and the Organizing Committee for the 
nice atmosphere created at this Symposium. Financial support from D.A.A.D. is 
gratefully acknowledged. 

7. References 

1. A. Bassetto, M. Ciafaloni, G. Marchesini, Phys. Rep. 100 (1983) 202 

2. Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller, S.I. Troyan, Rev. Mod. Phys. 60 
(1988) 373; Basics of Perturbative QCD", Editions Frontieres, Gif-sur-Yvette 
CEDEX-France (1991) 

3. Yu. L. Dokshitzer, V. S. Fadin, V. A. Khoze, Phys. Lett. B115 (1982) 242; 
Z. Phys. C15 (1982) 325; A. Bassetto, M. Ciafaloni, G. Marchesini, A. H. 
Mueller, Nucl. Phys. B207 (1982) 189; Ya. I. Azimov, Yu. L. Dokshitzer, 
V.A. Khoze, S.I. Troyan, Z. Phys. C27 (1985) 65; Z. Phys. C31 (1986) 213 

4. TASSO Coll., W. Braunschweig et al., Z. Phys. C47 (1990) 187 

5. OPAL Coll., M. Z. Akrawy et al., Phys. Lett. B247 (1990) 617 

6. L3 Coll., B. Adeva et al., Phys. Lett. B259 (1991) 199 

7. ALEPH Coll., D. Buskuhc et al., Z. Phys. C55 (1992) 209 

8. DELPHI Coll., P. Abreu et al., Phys. Lett. B347 (1995) 447 

9. ZEUS Coll., M. Derrick et al., Z. Phys. C67 (1995) 93; 

10. HI Coll., S. Aid et al., Nucl. Phys. B445 (1995) 3; 

11. S. Lupia and W. Ochs, preprint MPI-PhT/95-86, September 1995, |hep- 
ph/9509249| ; more details in paper eps0803 submitted to the "EPS-HEP 95" , 
Brussels, 27.7-2.8.1995 

12. C. P. Fong and B. R. Webber, Nucl. Phys. B355 (1991) 54 

13. Yu. L. Dokshitzer, V. A. Khoze and S. I. Troyan, Int. J. Mod. Phys. A7 
(1992) 1875 

14. DASP Coll., R. Brandelik et al., Phys. Lett. B67 (1977) 358 

15. ADONE Coll., B. Esposito et al, Lett. Nuovo Cimento 19 (1977) 21 

16. MARK I Coll., J. L. Siegrist et al, Phys. Rev. D26 (1982) 969 

17. Yu. L. Dokshitzer, D. V. Shirkov, Z. Phys. C67 (1995) 449 

18. DELPHI Coll., CERN preprint CERN-PPE/95-28, paper eps0539 submitted 
to the "EPS-HEP 95", Brussels, 27.7-2.8.1995 



11 



